Luttinger liquid of trimers in Fermi gases with unequal masses 
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We investigate one dimensional attractive Fermi gases in spin-dependent optical lattices. We show 
that three-body bound states - "trimers" - exist as soon as the two tunneling rates are different. We 
calculate the binding energy and the effective mass of a single trimer. We then show numerically 
that for finite and commensurate densities rif = n±/2 an energy gap appears, implying that the 
gas is a one-component Luttinger liquid of trimers with suppressed superfluid correlations. The 
boundaries of this novel phase are given. We discuss experimental situations to test our predictions. 
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Recent advances with ultra-cold atoms are opening 
new prospects to address fundamental theoretical issues 
in direct experiments A long standing problem is 
whether superconductivity can coexist with the presence 
of an unequal number of up and down fermions. An in- 
triguing possibility is the celebrated Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state 0, where the superconduct- 
ing order parameter becomes modulated in space. The 
experimental search for polarized superfluids in atomic 
quantum gases has so far been restricted to 3D configu- 
rations [11,13] ■ A new and promising direction is to confine 
atoms in higly elongated traps, where the FFLO state is 
known to be very robust [5j-l7|, as confirmed by detailed 
numerical simulations Q. 

Another exciting topic that is currently being explored 
experimentally is the pairing in Fermi gases with unequal 
masses, like mixtures of 6 Li and 40 K near a heteronu- 
clear Feshbach resonance [jl-Hlj]. Alternatively, one can 
also trap a two-component Fermi gas in a spin depen- 
dent optical lattice so that the corresponding effective 
masses are different [l2j]. Assuming that the transverse 
motion of atoms is frozen by a strong radial confinement, 
the system can then be described by the ID asymmetric 



Fermi-Hubbard 13-15] model 



H. 



aH 



where U < is the on-site attraction and t a are the 
spin-dependent tunneling rates. Here Ci a annihilates a 
fermion with spin a at site i and h ia is the local density. 

For tf — t± Eq. (p} reduces to the Hubbard model. 
For £|_ = Eq. ([T]) is nothing but a spinless version of 
the Falicov-Kimball model [16[ , originally devised for ex- 
plaining metal-insulator transition in mixed- valence ma- 
terials via hybridization of itinerant electrons (f) and 
localized holes (4-). There is a devil's staircase struc- 



ture in the ground state of the system [17|. For finite 
ti Eq. ([T]) is often referred to as an extended Falicov- 
Kimball model, which has been recently suggested to 
feature Bose-Einstein condensation of d-f excitons with 
macroscopic polarization [l8| and a spontaneous ferro- 
electric state fl9j | . 



For equal masses, t\ = fcj., the exact (Bethe ansatz) 
solution shows that n-body bound states with n > 2 are 
generally forbidden 20]. When the tunneling rates are 
different, t± < f-f, the above model is no longer intc- 
grable and many interesting questions arise: are three- 
body bound states (trimers) allowed ? What are their 
properties ? Can trimers open an energy gap like pairs 
do ? And if so, what are the differences between the two 
gapped phases ? The purpose of this Letter is to pro- 
vide an explicit answer to these relevant questions. We 
show that there is formation of three-body bound states 
of two heavy (4,) fermions and one light (f) fermion. We 
first find the regime where these trimers are stable as a 
function of mass asymmetry and attraction strength by 
solving the three-body problem. We then use the DMRG 
method to show that at low but finite density there ex- 
ists a novel phase with a nonzero energy gap which is 
a one-component Luttinger liquid of trimers, with ex- 
ponentially suppressed superconducting FFLO correla- 
tions. Finally, we calculate the boundaries of the trimer 
phase in the grand canonical phase diagram. These re- 
sults are in agreement with the generic bosonization anal- 
ysis of Ref. |2l| where the role of higher harmonics of the 
density operator was elucidated. The trimers discussed 
here are a cold atom an alog of the trions recently ob- 
served in semiconductors [221 ] - Three-body bound states 
(though of different origin) have also been predicted to 
occur in three-component Fermi gases [23[. In the fol- 
lowing we fix the energy scale in Eq. (JTJ) by setting tf=l 
and assume tj, ^ 1 without loss of generality. 

Three-body problem. — We start by calculating the bind- 
ing energy and the effective mass of the trimer. The 
Schrodinger equation can be conveniently rewritten in 
integral form by using Green functions. For the three- 
body problem in a lattice one finds [241 ] : 
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where £(k, q) = ej.(fc) + ej.(g) +£f(-P — k — q), P being the 
quasi-momentum of the trimer and e CT (fc) = 2£ CT (1 — cos k) 
the energy dispersions of the two components. Moreover 



2 




FIG. 1: Solid lines: Binding energy B tf of a trimer (one light 
and two heavy fermions) as a function of the hopping ratio t± 
for increasing values of the attraction strength U = — 2 (bot- 
tom), — 4, — 8, — oo calculated from the exact solution of the 
3-body problem. Open circles represent DMRG calculations 
in a chain of L = 100 sites. Inset: Inverse effective mass of the 
trimer as a function of the hopping ratio tj. for increasing val- 
ues of the attraction strength U = —2 (top), —4, —8, shown 
by solid lines. At the symmetric point t± — 1 the bound state 
disappears and the effective mass reduces to the sum of the 
masses of constituents (dashed lines). 



R E (q) = [l + UI E {q)} 1/2 , with 
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Eq. |2]) can be considered as an eigenvalue problem 
Kg • f = Af , where the energy E is fixed by the constraint 
A = 1. We solve this equation numerically for zero quasi- 
momentum P = 0. The binding energy E\ x of the trimer 
is related to the total energy E by —E = £p air + E\ T , 
where £^ air = -2(1 + t i ) + ^U 2 + 4(1 + t^) 2 is the bind- 
ing energy of the constituent pair [25[. In Fig. [T]we plot 
the binding energy of the trimer as a function of the mass 
asymmetry t± for increasing values of the attraction U. 
We see that E^. vanishes at the symmetric point tj, = 1, 
in agreement with the Bethe Ansatz solution 20|. As 



the mass asymmetry increases the binding energy also in- 
creases until it saturates at ij, — > 0. In this limit the func- 
tion reduces to a constant ifi(fc) = 1/ y/E{E — 4), im- 
plying that Eq. ([2]) has solution of the form f(k) = sinfc. 
Substituting this into Eq. ^ and integrating over mo- 
mentum we find E = — U 2 /(l — U). This gives 



El(h = o) = 



u 2 



1 - u 



+ 2 - yJU 2 +4 , 



(4) 



in agreement with our numerical solution. In particular, 
for infinitely strong attraction, Eq. ((4]) yields E\ x = 1, 
showing that the trimer binding energy remains Unite 
in contrast with the pair binding energy which is in- 
stead divergent. Indeed, when the heavy particles are 
at neighboring sites, the light fermion can hop from one 



site to the other without changing the interaction en- 
ergy. Therefore the total energy gain is at most equal 
to tf. In the strong coupling limit, \U\ ~> 1, Eq. @ 
can be solved by the ansatz f(k) = sink/R E (k) yielding 
E^ r (U = — oo) = (tj, — l) 2 , which is shown in Fig. Q] with 
black line. 

Let us now briefly discuss the effective mass M t * of the 
trimer which is defined by 1/M t * = d 2 E/dP 2 evaluated 
at P = 0. The inverse effective mass is plotted in the 
inset of Fig. Q] as a function of the hopping rate t\, and 
for different values of the attraction strength. We see 
that the trimer becomes heavier as t\, decreases or \U\ 
increases. At the breaking point, t± = 1, the effective 
mass reduces to (\/Mt.i. + l) 2 + U 2 + 2)/4t\., correspond- 
ing to the sum of the masses of a pair [25( , and of a heavy 
fermion. This quantity is plotted in the inset of Fig. Q] 
with dashed lines. 

Trimer gap. — We now turn our attention to the effects 
of trimers at finite density. We first calculate the trimer 
gap, namely the energy needed to break a single trimer. 
This is defined as 

A tr = - lim [E L (Ni + l,N± + 2) + E L (N t , 

L— >oo 

-£ i (JV t + l,2V; + l)- J B i (JV t ,JV4. + l)], (5) 

where El(N^, N±) is the ground state energy of a gas 
with spin populations N-f.Ni in a chain of size L. The 
limit in Eq. (|5|) is taken assuming Nu — > oo with n a = 
N a /L being fixed. The trimer gap is the generaliza- 
tion of the binding energy at finite density, with the ref- 
erence state being the many-body state (Nf,N±) rather 
than the vacuum (0,0). We evaluate Eq. (O numeri- 
cally via DMRG technique on lattices of up to L = 160 
sites with open boundary conditions and perform careful 
finite-size scaling in order to extract the thermodynamic 
limit behavior. For equal masses, t± = 1, our results 
are consistent with A tr = for any concentration, as 
expected from the Bethe Ansatz solution. For unequal 
masses, corresponding to =/= 1, the trimer gap ([5]) is 
finite only when the two concentrations are commensu- 
rate, namely = 2n^. 

Fig.[2]shows typical results for A tr at fixed t\, = 0.3 and 
n i/ n -f = 2. For vanishing density the trimer gap reduces 
to the binding energy E^ r , shown in Fig. [2] by arrows. As 
the density increases, A tr decreases and eventually van- 



ishes at a critical concentration 



in a sharp 



contrast with the case of equal densities, rif = n^, where 
the associated pairing gap is always positive for any fill- 
ing. In other words, the opening of the trimer gap is a 
non perturbative effect requiring finite coupling strength, 
or equivalently, low enough densities. This is consistent 
with our bosonization approach (2lj . 
Correlation functions. — The opening of the trimer gap 
drastically affects the ground state properties of the 
gas, since correlation functions of all operators breaking 
trimers fall off exponentially rather than algebraically. In 
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FIG. 2: Trimer energy gap plotted versus density n_j_ of the 
heavy component (n-f = n\/2) and for different values of 
the attraction U — —2 (bottom) and U = —4. The mass 
anisotropy is = 0.3. The gap Atr reduces to the bind- 
ing energy Et T at zero density (see arrows) and vanishes at 
a critical value of the density. The data are obtained from 
finite size scaling after DMRG simulations with system sizes 
L = 80, 120, 160, assuming a linear dependence in 1/L. In the 
inset we show the scaling analysis for two different concentra- 
tions nj. = 0.45 and n_i = 0.6 with U = —4. 



particular, the superconducting correlations decay as 



exp(-|* 



co S (Q|i-i|) (6) 



= \kl 



is 



where the decay length £ oc A^. , Q 
the FFLO momentum, and a is a nonuniversal num- 
ber. Two-point correlations (cl^Cja) display similar be- 
havior. In Fig. |3] we explicitly check this prediction by 
showing the superconducting correlations in the gapped 
(ni < nj r ) and gapless (n^ > nj r ) phases. We see that in 
the gapped phase the typical FFLO modulation is pre- 
served, as discussed in Ref. [13, EH , but quasi long range 
order is lost. 

Phase diagram. — The presence of trimers and other 
bound states induced by the mass asymmetry changes 
the topology of the grand canonical phase diagram of 
the gas. The latter is obtained by replacing the den- 
sities rif i by two new variables, corresponding to the 
mean chemical potential /i = dE/d(N^ + Ni) and the 
effective magnetic field h = dE/d(N^ — N\), where E 
is the ground state energy. The evolution of the overall 
shape of the phase diagram with changing has been 
presented in Ref. [la ]. Here we concentrate on the non- 
trivial changes due to the extra bound states induced 
by the mass asymmetry, which were not discussed previ- 
ously. To that end we work at fixed values of U and tj,. 
We approximate the derivatives by finite-difference for- 
mulas similar to Eq. (J5]) , and obtain the phase diagram 
shown in Fig. [4j The partially polarized phase (PP) cor- 
responds to configurations with imbalanced spin popula- 
tions (1 > nj, > rif > 0). This phase is limited from the 
right by the phase of equal densities (ED) (n^ = n-\-) and 
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FIG. 3: Superconducting correlations as a function of the 
distance from the center of the chain. The upper curve cor- 
responds to nj_ = 0.7, where the trimer gap is zero, cf. Fig. 
[2] The lower curve refers to low density nj. = 0.3, where 
A tr > 0. The parameters used are U = —4, t± — 0.3 and 
L = 200 and the densities are commensurate, = nj_/2. 
Notice the change from algebraic decay (upper curve) to ex- 
ponential decay (lower curve). 



from below by the fully polarized (FP) phase where the 
minority component is absent (n^ = 0) [2r3 |. For equal 
masses, ij. = 1 the three phases meet at a single point. 
In presence of n-body bound states with n > 2 this spe- 
cial point splits into an extended line, corresponding to 
a direct boundary between PP phase and vacuum. To 
see this, suppose there exists a bound state made of p t- 
fcrmions and q 4,-fermions, where p and q are non negative 
integers. Since the system density is zero, the Taylor ex- 
pansion El =oc (p, q) = (p + q)/i + (p — q)h becomes exact, 
yielding a straight line in the (h, /i) plane associated to 
the bound state. The true phase boundary fj, = Hvacih) 
with the vacuum is given by 



E L= 00 (p,q) -{p-q)h 

p + q 



(7) 



resulting in a piecewise straight line, cf. Fig. 2J While 
for equal masses the only states entering ([7]) are a sin- 
gle ^-fermion and a pair (p = 0, q = 1 and p = q = 1, 
respectively), for tj_ ^ 1 additional bound states appear, 
e.g. trimers (p = 1, q = 2), quadrimers (p — 1, q = 3), 
pentamers (p — 2, q — 3), etc. — leading to an extended 
PP-vacuum boundary, as shown in Fig. 0] It is also in- 
structive to consider the locus of n-j- = %/2 on the phase 
diagram. At low density (n^ < u£ r ), the trimer gap is 
non-zero and the commensurate densities occupy a finite 
area of the (h, /i) plane. At higher density (n^ > n"'), the 
energy gap closes and the locus shrinks to a single line, 
as illustrated in the inset of Fig. [4] As the mass asym- 
metry increases, the trimer phase grows in size. Similar 
behavior is found for the gapped phases associated to all 
other bound states (not shown in Fig. [4]). We emphasize 
that in all the simulations reported above we observed a 
uniform ground state, apart from the usual Friedel oscil- 
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lations due to the open boundary conditions. Collapsed 
phases occur for larger values of the ratio |E/|/i|, [14] . 

Relevance for experiments. — Let us finally discuss how 
our predictions can be tested in experiments with 
trapped Fermi gases. The asymmetric Hubbard model 
(UJ can be directly implemented by using spin-dependent 
optical lattices [LJj. The presence of shallow harmonic 
traps V„{z) = niLUz a .z 2 /2 can be taken into account via 
local density approximation, starting from the homoge- 
neous solution. Here m is the atom mass and U) Z(T the 
trapping frequencies. In order to form trimers, we require 
T, h\o za < E\ x , T being the temperature of the gas. For 
instance, consider a sample of 6 Li atoms in a lattice with 
periodicity d — 250 nm and tunnelling rates t-j- = 345nK 
and ti = 0.2tf. Assuming U = — 4if, from FigUJ we 
obtain E\ T = 0.3^ = 104 nK, well within experimental 
reach. The binding energy of trimers ( Fig . [1]) can be di- 
rectly measured by rf spectroscopy [27|, l28[ . The effective 
mass (inset in FigfT]) can be obtained by measuring the 
frequency of the dipole oscillations of the cloud [l|. For 
small displacements around the equilibrium position, the 
kinetic energy of trimers is quadratic, P 2 /2Af t * , and the 

latter is simply given by Wdip = y(^| + 2w|)m/M t *. To 
enter the trimer phase at finite density, the above con- 
dition becomes more stringent, namely T,hw ZG < A tr . 
Moreover both components must be degenerate, imply- 
ing T < Ep-fjEpi, where Ep a are the Fermi energies in 
the absence of interaction. This sets a lower bound on 
the values of the densities at the center of the trap. For 
the above choice of parameters, the best trade-off occurs 
around 7^(0) = 2ra t (0) - 0.3/d yielding T,huj z<7 < 60 
nK. Notice that to have an extended trimer phase, the 
two spin populations should be tuned close to the com- 
mensurate point Nf = Ni/2, otherwise phase separation 
in shells will occur 0] ■ Finally, the suppression of the 
superconducting correlations (FigfS]), signalling the emer- 
gence of the trimer phase, can in principle be detected 
using interferometric techniques, as discussed in Ref . [29jj . 

Concluding, we have shown that ID attractive fermions 
with unequal masses form trimers and other more exotic 
bound states. Differently from pairs, these states can 
only open a gap at low density or, equivalently, strong 
interactions. In the gapped phase FFLO superconduct- 
ing correlations are exponentially suppressed. The prop- 
erties of trimers in vacuum and at finite densities are 
experimentally accessible with ultra-cold atoms. Numer- 
ical simulations were performed using the ALPS libraries 
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